New Strong-Field QED Effects at ELI: Nonperturbative Vacuum Pair Production 
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Since the work of Sauter, and Heisenberg, Euler and Kockel, it has been understood that vac- 
uum polarization effects in quantum electrodynamics (QED) predict remarkable new phenomena 
such as light-light scattering and pair production from vacuum. However, these fundamental effects 
are difficult to probe experimentally because they are very weak, and they are difficult to analyze 
theoretically because they are highly nonlinear and/or nonperturbative. The Extreme Light Infras- 
tructure (ELI) project offers the possibility of a new window into this largely unexplored world. 
I review these ideas, along with some new results, explaining why quantum field theorists are so 
interested in this rapidly developing field of laser science. I concentrate on the theoretical tools that 
have been developed to analyze nonperturbative vacuum pair production. 



I. INTRODUCTION 

The Extreme Light Infrastructure (ELI) project will provide lasers with electromagnetic fields of unprecedented 
intensity, allowing a new experimental window into the largely unexplored regime of nonperturbative quantum elec- 
trodynamics (QED). This has implications not just for QED, but also for fundamental issues in quantum field theory, 
as well as nuclear, atomic, plasma, gravitational and astro- physics. Quantum vacuum fluctuations mean that the 

@ED vacuum behaves like a polarizable medium that modifies classical behavior, leading to novel quantum effects 
, H i, i, H . Many of these, such as the Lamb shift BiBl, Delbriick scattering [§| , photon splitting [lo|, [lH , nonlinear 
Compton scattering [l^ , and the Casimir effect [isl . ll^Tlsj , have been experimentally observed, while others such as 
elastic photon-photon scattering and nonperturbative electron-positron pair production from vacuum, have not yet 
been observed. This is essentially because they are very weak effects. Our current understanding of the perturbative 
regime of QED is extremely good. The classic example is the anomalous magnetic moment g of the electron [defined 
by the proportionality between magnetic moment and spin: fl — g^S], concerning which there has been dramatic 
recent progress, both experimentally and theoretically. On the theoretical side a four-loop computation [involving 
891 four-loop Feynman diagrams!] gives g as a series in the fine structure constant a — y- [which sets the interaction 
strength of perturbative QED] as: 
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On the experimental side (l7j . a single-electron cyclotron has been used to measure g directly, without measuring a, 
and the result agrees remarkably well with the value obtained from (jl.ip using a precise independent measurement of 
a (from atomic recoil experiments with Rb): 
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This is an impressive confirmation of the precision of perturbative QED. 

By contrast, we know very little about the nonperturbative regime of QED that arises when we consider QED in 
ultra-strong external fields. To quantify "ultra-strong", recall the computation of Heisenberg and Euler [3] of the 
probability of vacuum pair production in an applied uniform electric field of strength £. Vacuum pair production is 
a process in which virtual dipole pairs in the vacuum can be accelerated apart by the external field, becoming real 
asymptotic e"*" e~ pairs if they gain the binding energy of 2mc^ from the external field, as depicted in Figure [TJ This 
is a non-perturbative process, and the leading exponential part of the probability, assuming a constant electric field, 
was computed by Heisenberg and Euler 0, 0] : 
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FIG. 1: Pair production as the separation of a virtual vacuum dipole pair under the influence of an external electric field. 



building on earlier work of Sauter [l8|. This result sets a basic scale of a critical field strength and intensity near 
which we expect to observe such nonperturbative effects: 
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As a useful guiding analogy, recall Oppenheimer's computation ^19i] of the probability of ionization of an atom of 
binding energy E'b in such a uniform electric field: 
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Taking as a representative atomic energy scale the binding energy of hydrogen, _Bf, 
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13.6 eV, we find 



(1.6) 



This result sets a basic scale of field strength and intensity near which we expect to observe such nonperturbative 
ionization effects in atomic systems: 
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These, indeed, are the familiar scales of atomic ionization experiments. Note that ^ionization (jjffgj-g from £c by a factor 
of ~ 4 X 10~^. These simple estimates explain why vacuum pair production has not yet been observed - it is an 
astonishingly weak effect with conventional lasers [201 . This is because it is primarily a non-perturbative effect, 
that depends exponentially on the (inverse) electric field strength, and there is a factor of ~ 10'' difference between 
the critical field scales in the atomic regime and in the vacuum pair production regime. Thus, with standard lasers 
that can routinely probe ionization, there is no hope to see vacuum pair production. However, recent technological 
advances in laser science, and also in theoretical refinements of the Heisenberg-Euler computation, suggest that lasers 
such as those planned for ELI may be able to reach this elusive nonperturbative regime. This has the potential to open 
up an entirely new domain of experiments, with the prospect of fundamental discoveries and practical applications, 
as are described in many talks in this conference. 



II. THE QED EFFECTIVE ACTION 

In quantum field theory, the key object that encodes vacuum polarization corrections to classical Maxwell electro- 
dynamics is the "effective action" T[A\, which is a functional of the applied classical gauge field A^{x) [2^ [23l. |2^. 
The effective action is the relativistic quantum field theory analogue of the grand potential of statistical physics, in 
the sense that it contains a wealth of information about the quantum system: here, the nonlinear properties of the 
quantum vacuum. For example, the polarization tensor II^i^ = ^ contains the electric permittivity and the 
magnetic permeability /iy of the quantum vacuum, and is obtained by varying the effective action r[^] with respect 
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to the external probe Af^{x). The general formalism for the QED effective action was developed in a series of papers 
by Schwinger in the 1950's [12, [2^. T[A] is defined [23j in terms of the vacuum- vacuum persistence amplitude 
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(2.1) 



Note that T[A] has a real part that describes dispersive effects such as vacuum birefringence, and an imaginary part 
that describes absorptive effects, such as vacuum pair production. Dispersive effects are discussed in detail in Holger 
Gies's contribution to this Volume. The imaginary part encodes the probability of vacuum pair production as 
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Here, in the last [approximate] step we use the fact that Im(r)/?i is typically very small. The expression (12. 2p can 
be viewed as the relativistic quantum field theoretic analogue of the well-known quantum mechanical fact that the 
ionization probability is determined by the imaginary part of the energy of an atomic electron in an applied electric 
field. 

From a computational perspective, the effective action is defined as [l^, H^l 



r[^] h In det [ip- m] 
— ?itr In [ip— m] 
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Here, p = j^^Df^, with 7^ being the Dirac gamma matrices of relativistic quantum mechanics, and the covariant 



derivative operator. Da ~ du 



-A^, defines the coupling between electrons and the electromagnetic field A^. When 



the gauge field A^ is such that the field strength, = d^j^A^, — d,jA^, is constant, the spectral problem becomes that 



of 2 harmonic oscillators, as is easily seen in the Fock-Schwinger gauge: 



In this case, F^^^, is a 4 x 4 



antisymmetric matrix, so its eigenvalues come in ± pairs, and can be expressed simply in terms of the two relativistic 
invariants {E^ — B^) and E ■ B. Then the determinant can be computed in closed form 0, Q. This is like a relativistic 
analogue of Landau's famous computation of diamagnetism for nonrelativistic electrons in a constant magnetic field 
PH . For a constant electric field, of strength the result of Heisenberg and Euler is 
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The poles of cot(s) lead to an imaginary part of F, whose leading behavior is determined by the first pole at s = tt: 
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(we now work in units with fi — c— 1). This result, and its interpretation in terms of vacuum instability, was stated 
clearly already by Heisenberg and Euler building on earlier ideas of Sauter Il 8|. 

There is an interesting interpretation of this result as a virial expansion [26l . l27l l28j , in which we write 



-ImF = - ^trp ln(l - 



(2.6) 



where Up = exp [— 7g-("^^ ^^e mean number of pairs with momentum p±^ transverse to the direction of the 

field. The sum is over spin states s, and the momentum phase space trace becomes: trp = Vj_ J (^Jj^^-i 1 which 

we recognize the electric field analogue of the familiar Landau-level degeneracy factor. Typically, Up is very small, 
so we can approximate — ln(l — np) « rip, in which case we immediately obtain the leading Heisenberg-Euler result 
(|2.5p . On the other hand, expanding the logarithm we obtain the full instanton sum of Schwinger [23.] 
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This instanton sum expression can also be obtained by summing over all the poles at s = n7r in (|2.4p . 
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This result is for spinor QED. If we were to consider the quantized charged particle to be a scalar particle, or even 
to have classical statistics, then the analogous expressions would be: 

{~ Y.S t^'p ~ ^*p) ' Fermi 
trpln(l + rip) , Bose (2.8) 
trp Up , Boltzmann 

Thus, when Up is small, the leading term is independent of statistics [apart from a simple degeneracy factor (2s + 1)], 
and so many computations work with scalar QED rather than spinor QED, since it is notationally simpler. But we 
should remember that ^-ImT « trp rip is an approximation. The higher order terms in the sum (|2.7p correspond to 
the coherent production of multiple pairs [26i, i27i. i28] in the same spacetime volume A^, and quantum statistics effects 
become important. It would be very interesting if these higher order correlations could be probed experimentally, but 
at present we are concentrating on being able to measure directly the first term. 



III. BEYOND THE CONSTANT FIELD APPROXIMATION: ONE-DIMENSIONAL 

INHOMOGENEITIES 



So far, we have discussed the theoretical computation of the vacuum pair production probability only in the 
approximation of a constant and uniform applied electric field. But to reach the nonperturbative regime of £ ^ £c, 
for example with ELI, will involve intense, short-pulse, focussed laser beams. Therefore, we need to understand how 
the Heisenberg-Euler result (|2.5p for Im(r[A]) is modified for gauge fields Af^{x) corresponding to reahstic laser fields 
with strong spatial and temporal inhomogeneities, as depicted in Figure [2l This is a well-posed problem, for which 
the required formalism exists [2^ . [23l . [2^ . in terms of determinants of the Dirac operator. It is nevertheless a very 
difficult computational problem for which there are few quantitative results. A somewhat analogous situation exists 
in the theory of the Casimir effect, as sketched in Figure [31 Casimir's original computation jl^ for a configuration of 
two parallel plane perfect mirrors leads to a simple expression for the Casimir energy and force, but for more realistic 
experimentally realizable surface configurations the computation is well-defined and easily posed, but it is extremely 
difficult to obtain explicit quantitative results [isj . 




FIG. 2: Progression from Heisenberg and Euler's original configuration of a uniform constant electric field, to more realistic 
configurations with one- dimensional inhomogeneities, and eventually to realistic focussed laser pulses. Here, non-perturbative 
aspects of vacuum polarization are probed by the externally imposed electromagnetic fields. 



o 




FIG. 3: Progression from Casimir's original configuration of perfect parallel plates, to realistic experimental configurations with 
difi'erent geometries and imperfect surfaces. Here vacuum polarization is probed by the externally imposed boundaries. 
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A. The Keldysh approach 



As a first step in this direction of computing the vacuum pair production probability in inhoniogeneous external 
fields, recall once again the analogy between vacuum pair production and atomic ionization. In the context of atomic 
ionization, Keldysh (soj considered a monochromatic time dependent electric field represented by the sinusoidal func- 
tion £{t) = £ cos{u!t). For an atom with binding energy iSf,, the peak electric field amplitude £ defines a frequency scale 
iVfc = e£/y/2mEb, and Keldysh identified two important physical regimes characterized by the "Keldysh adiabaticity 
parameter" 



1K= — = ^ ■ (3.1) 

luk e£ 

The region ^ 1 is an instantaneous, tunneling, nonperturbative regime, while the opposite limit where ^ 1 
turns out to be a perturbative multi-photon regime. Note that both regions, small and large 7^, can be spanned 
without violating the relevant approximations: 
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£ < £i°ni''atio„ ^ V j^^^j^ g^j^ approximation] 

en 



huj <^ Eh [classical field approximation] (3.2) 

In the sinusoidally time dependent electric field, £{t) = £cos{ujt), Keldysh found the following simple approximate 
expression for the ionization probability: 



-^ionization ^ exp 
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The function .9(7!-) ^ 1 for jk ^ 1, while (7(7^) ^ 3\n{2'-fK) / {2jk) for •jk ^ 1- Thus, we obtain the limiting 
behaviors: 

/' r /2 " 

, "fK ^ 1 (nonperturbative) 
, ^ 1 (perturbative) 

For small 'Jk the field £{t) is essentially constant over the time scale I/uik, and the ionization process is a nonper- 
turbative tunneling effect, as in Oppenheimer's result (jl.Sp . On the other hand, for large 7^ it is a perturbative 
multi-photon effect, which can be seen from the fact that in this case the probability is given by the normalized 
perturbative field £ raised to a power equal to (twice) the number of photons required to reach the binding energy. 
Keldysh's work was developed further by several authors [l^, [SJ] , and the resulting theory to describe strong- 
field ionization in atomic and molecular systems is commonly referred to as the "ADK" theory (see also 36]). 
However, an important comment should be kep t in mind: in ultra-intense lasers it is not necessarily realistic to neglect 
the magnetic component of the laser field [33], so that the representation of the laser field by just a time-dependent 
electric field £{t) should be extended. We return to this comment below, in Section [TVl 

The Keldysh approach has been adapted to the QED vacuum pair production problem [ssl . [39l . [40| , once again with 
the laser field modeled by an oscillating electric field £{t) = £cos{ujt), and with the role of the binding energy Et, 
played by 2mc^. The frequency scale set by the peak electric field amplitude is luqed — e£/{mc), and the analogous 
" adiabaticity parameter" is 

w mcuj 1 ,„ 
7= = — ^ = - . (3.4) 

UJQED eo ao 

Here we note that we can also express the adiabaticity parameter 7 as the inverse of the "normalized field strength" 
parameter 

ao = . (3.5) 

mcui 

The parameter ao provides a familiar characterization of the laser intensity in the laser and plasma physics literature. 
Using WKB techniques, one finds [sl,!!^!!^], analogous to 




^ ^ J , 7 < 1 (nonperturbative) 

-f^airprod. , ^.9rr,.^/t„. . ^'^■^> 

7 3> 1 (perturbative) 
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In the perturbative mu lti-p hoton regime, this QED pair production effect has been observed in a beautiful experi- 
ment (E-144) at SLAG |43|, in which a laser pulse collided with the (highly relativistic) SLAG electron beam, leading 
to nonlinear Gompton scattering involving 5 photons, producing a high energy gamma photon that decays into an 
electron-positron pair. An important challenge for the ELI project will be to observe pair production in the non- 
perturbative regime, and directly from vacuum, as a projected goal of ELI l'| is to attain an adiabaticity parameter 
7 ~ 1/5000 ^ 1. Theoretically, this requires reliable predictions of pair production rates and spectra in this non- 
perturbative regime, for specific ELI laser parameters. Before discussing new work in this direction I review various 
approaches that have been developed to analyze the case of an external field that is inhomogeneous in one direction. 



B. The one-dimensional scattering picture 



An important perspective on vacuum pair production in electric fields with a one-dimensional inhomogeneity is 
obtained by viewing the problem as a one dimensional quantum mechanical scattering problem (Popov [40| and 
Perelomov and Zel'dovich [s^l each attribute this idea in this context to L. Pitaevskii). This QM scattering problem 
can then be solved exactly in some cases, numerically in most cases, or semiclassically by various WKB approaches. 
Gonsider a time dependent electric field represented by the vector potential ^3(t), so that Esit) — ~As{t). As 
mentioned above, for the leading contribution to the pair production probability we can simplify matters by working 
with scalar QED, so we write the Klein-Gordon equation {—D^D^ + m'^)4> = for the field operator cj) as 

-4>-{P3-eA3{t))^=im^+pm (3.7) 

This can be viewed as a one dimensional quantum mechanical Schrodinger equation, with t playing the usual role 
of X. To describe the process of pair production, we seek solutions of ()3.7p with just positive frequencies in the far 
future, but a mixture of positive and negative frequencies as i — oo: 



bff ' ' 



FIG. 4: The Bogoliubov transformation between in and out states at t — ±oo, with boundary conditions in (|3.8|) . can be viewed 
as a one dimensional quantum mechanical scattering problem. 



~ ap-e-**V"'+p' , t^+oo (3.8) 

In the QM picture, this looks like a scattering problem, and the Bogoliubov transformation between t — ±oo modes 
is encoded in the reflection amplitude bp. Since the "potential", —jp^ — e^3(t))^, is negative, and the "energy", 
(m^ + p\), is positive, this is the case of over-the-barrier scattering [41 | . and so the reflection coefficient is 
expected to be exponentially small. By familiar manipulations (ssl . l39l |40|. |4^ , the pair production probability is 
given by 

-2 Im ^ y/m'^+pl + [P3-eA3{t)]^ dt 

(27r)3 

The dominant exponential contribution can be extracted by setting p = 0, while the prefactor comes from Gaussian 
integration in the momenta, after expanding the integral in the exponent to quadratic order in momenta. The 
transverse and longitudinal momenta contribute differently, and one arrives at a simple compact formula [39| : 



^5(7^) 
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(3.10) 
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Here the function 5(7^) is analogous to the function g{jK) 
formula for g is: 

ff(7^) = - / d,y 



appearing in Keldysh's expression (|3.3p . An explicit 



1/1 



(3.11) 



where we have written the imaginary time version of the gauge potential as ^ f{it), defined y = //7, and it is 
understood that /' is written back in terms of y. For example, for the cosine field £{t) = £cos{ujt), we have 
^3(0 ~ ~Zj sin(ajt), and so /(X4) = sinh(ci;a;4). Thus, with y = f /j, we have /' = cosh(wx4) = -^/l + ■ 



"f^y^, and so 



4^1 + 7^ 



1 - 
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1+7' 
7 < 1 
7 > 1 



E 
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1+7' 



(3.12) 



Here K and E are the standard complete elliptic integrals. The result in (|3.12p explains the tunneling and multiphoton 
limits quoted in (|3.6p . Notice that in the WKB formula (|3.10p . the exponent and the prefactor are both expressed in 
terms of the same function .9(7^). This result is typically obtained in the imaginary time method of WKB [39l. [40|. 
but can also be found in conventional WKB 'ss'l, and in the phase integral method [i^l, where higher order WKB 
contributions have also been considered. A similar formalism exists for static electric fields that depend on one spatial 
coordinate, represented for example by the gauge field ^0(^3) [i^, 111,14^. It is useful to compare these semiclassical 



results with some exact results [45|, |46| that arise for the cases £{t) = £ sech (ut) or £"(0:3) = fsech (kxs). For 
these fields, the associated Klein-Gordon equation is a soluble Schrodinger problem and so the effective action can be 
computed in closed form. (In fact, the associated Dirac equations are also exactly soluble). The semiclassical results 
compare very well with these exact results in the semiclassical regime (i^ . |47| . For example, using the phase integral 
method, for the soluble cases £{t) ~ f sech^(tji) or f (0:3) — £ sech^(fc 0:3), respectively, the particle number can be 
expressed as [il] 



M 



(1 + 7^) 



2\5/4 



exp 



(e£)^ 5/ 



exp 




where in the first expression 7 = while in the second expression 7 — When the gamma factors go to we 
recover the form of the constant field Heisenberg-Euler result, as is explained in detail in '44'|. These results illustrate an 
important point: a temporal inhomogeneity tends to enhance particle production, while a spatial inhomogeneity tends 
to suppress particle production. This follows from the different behaviors in the above formulas of the exponents and 
of the prefactors. Physically, this can be understood as follows: with a spatial electric field, the dominant nucleation 
of pair production occurs at a maximum of the field, but the virtual particles are then accelerated apart, away from 
this point, and if the field falls off too quickly then they may not gain enough energy to become real particles, so that 
the pair production is suppressed. This effect is seen clearly in the vanishing of the particle number when 7 — > 1. On 
the other hand, with a time-dependent field we can regard this as a tunneling problem with an oscillating barrier, for 
which the tunneling is enhanced roughly speaking because the average barrier is lower (33l . |40| . 



C. Quantum Kinetic Approach 



Another widely-studied approach in the analysis of the time dependent electric field case, £ — £{t), is the quantum 
kinetic approach [i^ . [49l [sol . [sil . [s^ . The key quantity in this formalism is the momentum distribution function, 
A/'(p, i), of pairs at a time t, which satisfies a quantum Vlasov equation (here for spinor QED): 

Here we have defined 



uj^{t) ^m^+pl + (p3 - eAsitjf 



(3.14) 
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For scalar QED the analogous quantum Vlasov equation is: 



dt 2uj{t) 



dt''^^(l + 2Afsc(t'))cos(2 [ dtLj(t)] (3.15) 
^{t') \ Jt' J 



The asymptotic large-time value of M{p, t) gives the momentum distribution of the produced pairs. Note that since 
the field is spatially uniform we can label the produced pairs by their momentum. The total particle number is 
obtained by integrating Af{p,oo) over the momenta. The form of these quantum kinetic equations shows that the 
process is inherently non-Markovian and non-local in time [i^. Furthermore, back-reaction effects have been studied 
in this approach |50|. showing that they are negligible if the maximum amplitude of the electric field is lower than 
the critical field £c [53| . 

It is not as widely appreciated that the quantum kinetic approach is equivalent to the Popov/Pitaevskii one 
dimensional scattering approach described in the previous section jsj]. We illustrate this briefly here for scalar QED. 
Recall that we can write the Klein-Gordon equation (|3.7p as 

^ + io^{t)(j) = (3.16) 

and convert this Schrodinger-like scattering problem into a Riccati equation as follows [40*1. First, express the complex 
field (j) as 

(j) = A {e-''^ -f Re'^P) ; if{t) = f uj{t')dt' (3.17) 



in terms of two other complex fields R and A where 

A Re^f 



(3.18) 



A e-^f + Re^f 
Then 

4> = -iujA{e-''f -re'^P) , (3.19) 
and the scattering equation (j3.7p becomes a Riccati equation for the scattering amplitude R: 

R=^ {e-^'"^ - R^e^'-^) . (3.20) 

It is straightforward to integrate this equation numerically, with initial condition R(—oo) = 0, and obtain the prob- 
ability for particle production as |i?(oo)p. Since (|3.7p has the form of a Schrodinger scattering problem with an 
inverted potential, this is the case of over-the-barrier scattering, which is conveniently analyzed with the method of 
imaginary time |55|. On the other hand, it could just as well be solved "exactly", by numerical integration. 

The connection with the quantum kinetic approach is seen by noting that in the quantum kinetic approach [for 
scalar QED], one writes the field in terms of two other complex fields a and (3: 

a 3 ■ 



^/2uj V2w 



e'^P . (3.21) 



^^/2uj \/2cJ 
The second equation requires a and (3 to be related by: 

2uj 2lj ^ ' 

Comparing the two decompositions for (j) in p.l7|) and (|3.21|) we see that R = [3 /a, and the equations for a and /? in 
p.22p are then equivalent to the Riccati equation p.20p . The particle number is 

A4c = = , (3.23) 

using the scalar relation: |ap — ~ 1. Converting the Riccati equation for the scattering amplitude R into 
an equation just for the probability one finds that it is in fact completely equivalent to the quantum Vlasov 

equation (|3.15p of scalar QED [54]. Thus, the quantum kinetic approach may be viewed as a real time version of 
Popov's imaginary time scattering problem. This can also be described equivalently in terms of Green's functions and 
the S-matrix [5^ , or in a Schrodinger picture [H^] . 
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D. Pulse shaping 

In going beyond the constant field approximation, the goal is to understand how inhomogeneities modify the constant 
field result (jl.3p . so that one might design a pulse shape that leads to a higher yield of vacuum pair production, at 
a lower peak field intensity. The importance of pulse shape for these processes has been studied widely |58i, ^ , and 
reflects experience gained in the strong-field atomic and molecular physics community, where sophisticated pulse- 
engineering techniques have been developed for ionization processes [60, [6l[ . 

In several recent papers [6^ . estimates based on taking time dependent laser pulse(s), neglecting the spatial variation 
in the focus or the interaction region, suggest that the intensity at which one could see the nonperturbative effects of 
vacuum pair production could be a factor of 100 smaller than the critical intensity /c « 4 x 10^^ W/cm^. This is once 
again reminiscent of the situation for atomic ionization, where ionization effects can be observed at intensities lower 
than the critical field estimate of lO^^W/cm^ coming from Oppenheimer's constant field result ()1.5p . In the QED 
vacuum polarization problem, with focussed laser pulses, the field polarization plays an important role (6^ . but the 
essential enhancement mechanism is that the interaction region is very large compared to the scale of the Compton 
wavelength, and so there is a significant prefactor enhancement. In these papers, the authors work with an external 
field that is an exact solution of Maxwell's equations, and that mimics a Gaussian laser beam; they then convolve 
this solution with a pulse time-envelope function to obtain an approximate solution that represents a focussed pulsed 
laser beam. In the vicinity of the focus spot, they argue that the spatial dependence can be neglected, so that in 
the end the problem becomes essentially one-dimensional, and the semiclassical results can be used. An important 
feature of this approach is that the role of the polarization of the field is clearly incorporated, and the cases of a single 
focussed beam, and of two counter-propagating beams were analyzed. The case of two counter-propagating circularly 
polarized beams leads to a higher yield of produced pairs, and the threshold intensity at which one might expect to 
see such effects is several orders of magnitude lower than the Schwinger limit. This is encouraging news for ELI. 

In a related development, it has recently been shown that the superposition of two time-dependent electric fields 
can lead to a significant enhancement of vacuum pair production, a " dynamically assisted Schwinger mechanism" [6^ . 
The specific example studied in [63] consists of a combination of a strong but slow field pulse, and a weak but fast 
field pulse: 

fsiow(i) = Ssech^ (fit) 
ffast(i) = esech^{ujt) 

The parameters fall into the hierarchies: 

< e < £ < £c 

Surprisingly, even though both field amplitudes, £ and e, are below the critical field £c in (|1.4p . there is significant 
enhancement of the pair production rate when the frequencies are also related as in (j3.24p . The non-perturbative pair 
production process that we would associate with the slow strong field interacts with the perturbative multiphoton 
pair production process that we would associate with the fast weaker field to produce a stronger impact than each 
process separately. This illustrates another possible approach to shaping the laser pulses in the time domain so as to 
enhance the pair production yield. 

IV. TOWARDS MORE GENERAL LASER FIELDS 

From a theoretical point of view, as well as from a practical experimental point of view, we should ask if we can go 
beyond the restriction to a one dimensional field inhomogeneity discussed in the previous section. There are several 
urgent reasons for studying more general background fields. First, such fields are needed to model more realistically 
the intense, focussed, short-pulse laser fields to be produced at ELI. Second, in the ultra-relativistic regime of ELI it 
is by no means clear that one can continue to neglect the magnetic component of the laser field and consider just the 
electric component. These questions need to be answered. 

In fact, surprisingly little is known about vacuum pair production in external fields that have a more complicated 
inhomogeneity structure, such as appears in a more realistic laser field. The results of all the QED work summarized 
in the previous section is that when the background field is approximated as an electric field pointing in one fixed 
direction, but with amplitude that either (i) varies with time, or [but not both] (ii) varies with space only along the 
direction of the field, then the problem can be reduced to a one-dimensional problem that can solved numerically, 
or approximately with WKB. However, all these standard approaches become practically intractable when the field 
depends on more than one coordinate. An interesting exception is for fields that depend on a light-cone coordinate 
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A. Worldline formalism for the effective action 



It is possible to develop semiclassical approximations for the case of multi-dimensional inhomogeneities in the 
external field, using Feynman's worldline formalism of the effective action. Feynman [66| formulated a first-quantized 
form of QED, which amounts to representing [63, HI] the effective action as a quantum mechanical path integral over 
closed loops Xfj,{T) in four dimensional spacetime, with the closed loops being parametrized by the proper time t. As 
Feynman noted, this proper-time representation is the natural representation for maintaining relativistic invariance. 
The propertime parametrization had been developed earlier by Fock and Nambu ^6S] , and was also used by Schwinger, 
in operator form rather than in path integral form, in his landmark QED computation of vacuum pair production 
[23| . Feynman's worldline path integral formalism has since been extended significantly, primarily for applications in 
perturbative quantum field theory ^68j , building on analogies and motivation from the Polyakov formulation of string 
theory. This rebirth has led to many beautiful advances in our understanding of perturbative scattering amplitudes 
[7l| . but here I propose to use it instead to extract non-perturbative information. Let us consider scalar QED rather 
than spinor QED, as it is notationally simpler, and the leading imaginary part of the effective action only differs from 
that of spinor QED by a spin degeneracy factor of 2. Then the effective action for a scalar charged particle (charge 
e, mass m) in a Euclidean classical gauge background A^(x) is the functional: 



r[A] = -tr ln(-i:>i^, 





°° dT 



T 







x(T)=x(0)=xW L ° \ 



(4.1) 



To understand the origin of this expression, recall the standard representation of the logarithm of an operator as: 
InM = — ^ e^^^ . This explains the second line in (|4.ip . where we have performed the trace in the position 

representation. The last line in (j4.ip follows from interpreting the matrix element Euclidean 
transition amplitude, which is then written using Feynman's representation as a path integral, but now involving paths 
x^ in four dimensional space. Since for the effective action we compute the trace of the logarithm, the matrix element 
is between identical points x^"^ , and so the path integral is over closed loops. Thus, we arrive at a first-quantized 
representation of the effective action r[A] in terms of a functional integral / Vx over all closed Euclidean spacetime 
paths x^{t) that are periodic (with period T) in the proper-time parameter r [66|. I use the QED worldhne path 
integral normalization conventions of (68j . Note also that we have expressed the path integral in Euclidean space, as 
is appropriate for describing a non-perturbative process. 

This representation of the effective action translates the primary computational problem into one of computing the 
quantum mechanical path integral in the last line of (|4.ip . This simple observation has several advantages. First, 
one could evaluate this Euclidean path integral numerically, using the worldline Monte Carlo technique of Gies and 
Langfeld [t^]. One creates a statistical ensemble of configurations, namely a set of closed loops Xfi{T) in 4 dimensional 
Euclidean spacetime, weighted by a factor exp[— drx^/A], and then one computes numerically, using a Monte Carlo 

algorithm, the expectation of the "Wilson loop" operator exp[— e J^^ drA^Xfj]. For the pair production problem, the 
exponentially small imaginary part, arising from poles in the T integral, can then be extracted numerically. This 
is still a challenging numerical problem, but it has been shown to work well for situations with one-dimensional 
inhomogeneities of the electric field |73j| . This is potentially a very powerful approach, as it does not rely on any 
particular symmetry of the background field - the ensemble of closed spacetime loops used in the numerical Monte 
Carlo evaluation of the quantum mechanical path integral is independent of the background field. Physically, the 
ensemble probes the structure and impact of the background field. This means that in principle the numerical worldline 
Monte Carlo approach can be used for quite general inhomogeneous laser field configurations. This approach deserves 
further exploration. 



B. Worldline Instantons 



A second, complementary, approach is more analytical: make a semiclassical approximation to the Euclidean 
worldline path integral. The essential idea is that the quantum mechanical path integral may be approximated by the 
contribution of a small number [possibly just one] of classical closed loops, plus the quantum fluctuations about these 
loops. These dominant Euclidean loops are called "worldline instantons". This is the Euclidean space analogue of 
the familiar semiclassical scattering theory (tS , [t^ , where a real-time scattering process is approximated by classical 
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trajectories, with corrections from quantum fluctuations about the classical trajectories. The dominance of the path 
integral by just a few closed paths is a semiclassical approximation, valid when the classical action is large compared 
to h. If the background field has a simple functional form, depending on just a few physical parameters such as peak 
field amplitude, a frequency or wavenumber, etc., then the regime of the semiclassical approximation can be precisely 
quantified in terms of a few dimensionless quantities. The idea of seeking dominant semiclassical solutions to the path 
integral appearing in (|4.ip was suggested long ago (76l |. and was first applied to the vacuum pair production problem 
for a constant electric field in More recently, the idea has been extended to the vacuum pair production problem 
with inhomogeneous background field configurations [78, 79, 80, 81]. The special dominant loops are solutions to 
the associated classical (Euclidean) equations of motion for a charged particle moving in the given [inhomogeneous] 
electromagnetic field Fi,„(x)\9(!^: 



(4.2) 



Solutions of the classical Euclidean equations of motion (|4.2p that are closed and periodic are called worldline instan- 
tons. 

For example, for a constant (Minkowski) electric field, £3=5, the Euler-Lagrange equations of motion in Minkowski 
space, with the coordinates x^{s) parametrized by proper time s are (here the overdot • refers to ^): 



2^3 



-2e£t 



-2 eE ±3 



The solution is well-known: 



cosh (2e£ s) 
2e£ 



t = ± 



sinh (2e£ s) 
2eE 



(4.3) 



(4.4) 



In a 1 + 1 dimensional Minkowski spacetime diagram, these follow open hyperbolic trajectories, as shown in Figure 
m On the other hand, in Euclidean space, with x^ = ict^ and s — s- ir, the equations read (here the overdot ■ refers 



X3 — —2e£x4 , X4 = 2e£ ±3 
and the solutions describe closed circular orbits, as shown in Figure [6j 



X3 ^± 



COS {2e£ t) 



X4 = ± 



sin {2e£ r) 



2e£ ' 2e£ 
For a time dependent (Minkowski) pulse field, given by £3{t) = £ sech^(ajt), the Minkowski equations are: 



(4.5) 



(4.6) 



FIG. 5: The Minkowski space-time trajectories of a charged particle in a constant electric field. The analogy with the trajectory 
of an accelerated observer is the starting point for the discussion of the Unruh effect in strong field laser systems - see R. 
Schiitzhold's article in this Volume. 



3:3 — ~2e£ sech^ {uj t) i 
t — ~2e£ sec\? {uj t) ±3 



(4.7) 
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FIG. 6: The Euclidean space-time trajectories [both directions around the loop] of a charged particle in a constant electric 
field. The open trajectories in Minkowski space, shown in Figure (5] are now closed trajectories in Euclidean space. This is the 
simplest worldline instanton, a periodic solution to the classical Euclidean equations of motion. 



The Minkowski solutions are the open trajectories: 

1 



1 



t{s) ~ ±— arcsinh 



arcsinh 



cosh(2ef Vl+7^s 



sinh 



(4.8) 



where 7 — These trajectories are shown in Figure [T] On the other hand, the Euchdean trajectories are closed 
loops: 



and are shown in Figure [3 



2^3 (t) = , 

2^4 (t) — — arcsin 



arcsinh 



7 cos (^2e£y/l + "f'^T 
, ^ sin (2e£\/l + 7^ r 



(4.9) 



FIG. 7: The Minkowski space-time trajectories [solid lines] of a charged particle in an electric field pulse, £{t) — £sec\)? {ujt), 
with adiabaticity parameter 7 = 2. 

Expanding in the fluctuations about the worldline instanton loop, x^{t) = a::J^'(r) +7y^(r), we obtain an approximate 
expression for the path integral: 



e / V-q exp 



2 7o 
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FIG. 8: The Euclidean space-time trajectories [both directions around the loop] of a charged particle in an electric field pulse, 
£{t) = f sech^(a't), with adiabaticity parameter ^ — 2. For reference, the dotted line shows the corresponding circular trajectory 
in a constant field where 7 = 0. The open trajectories in Minkowski space are now closed trajectories in Euclidean space. This 
is the simplest worldline instanton, a periodic solution to the classical Euclidean equations of motion. 



= , , (4.10) 

where S'[a;'^'] is the classical action, evaluated on the worldUne instanton loop. The prefactor is expressed in terms of 
the determinant of the fluctuation operator (the " Jacobi operator" ) is 

A^, = -^S^,-^+eF^4x^')^+edM^"')if ■ (4-11) 

In [tII it was shown that the worldline instanton exponential factor e^^^^°^^ agrees precisely with the WKB exponential 
factor exp[— ^^^j^ (7(7^)] in p.lOp . for the case of one-dimensional inhomogeneities. 

The form of the fluctuation operator in (|4.11|) exhibits another major advantage of the worldline approach - since 
the paths are parametrized by propertime, the fluctuation problem refers to a one-dimensional operator (actually, 
coupled one-dimensional operators) rather than a partial differential operator in four dimensional spacetime. This is 
technically important because there are simple mathematical methods for computing the determinant of such one- 
dimensional operators [zl, [H, [H, HI] , but these techniques do not generalize easily to partial differential operators. 
Thus one can numerically compute det(A) for any per iod T, and then perform the remaining T and x^^^ integrals 
in the effective action expression (14.11) . In fact, in (79| it was shown that for one-dimensional inhomogeneities the 
determinant of the fluctuation operator can be computed analytically: 

Det(A) = (2mifiT) T (4.12) 

V K(t)J / 

Further, after doing the T integral by steepest descents, one reproduces exactly the WKB result p.lOp . including the 
prefactor term [79|. 

This agreement with previous one dimensional semiclassical results is encouraging, but the real advantage of the 
worldline instanton method is that it can be extended to more general fields. Technically, there are several steps to the 
computation: (i) find periodic solutions [the worldline instantons] to the Euclidean classical equations of motion (|4.2p 

; (ii) The dominant exponential factor in the pair production rate is e~^^^ \ involving the classical action evaluated 
on the worldline instanton path; (iii) The prefactor coming from the semiclassical approximation to the quantum 
mechanical path integral is given by the determinant of the fluctuation operator in (|4.1ip , evaluated on the worldline 
instanton path; (iv) The T integration flxes the conserved quantity jif^ to take the value m^, and also produces 
a prefactor; (v) Possible residual dependence on the spacetime location of the loop must be integrated over. In a 
Gaussian approximation this also leads to another prefactor. 
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C. Worldline Instantons and the Gutzwiller Trace Formula 



Note that in general there are three apparently separate [but below we shall see that they are in fact related!] 
prefactor contributions to the final answer: one comin g fro m each integration in (j4.ip . This procedure is reminiscent 
of that used to derive the Gutzwiller trace formula 

[US [13 

for the trace of the Green's function in non-relativistic 

quantum mechanics: 



tvG{E) = J (fix dt 



^* (x\e^"'\x) 



„ e-kSp{E)-i^ nip 



(4.13) 
(4.14) 



Here Tp is the period of the p* classical closed orbit of energy E, Sp{E) is its action, Jp is the associated monodromy 
matrix [see below], and are Maslov indices. Before proceeding to the corresponding worldline instanton expression, 
we briefly sketch the strategy of the derivation of the Gutzwiller trace formula [an excellent introduction for physicists 
is in [ST*]]. The first step is to make a semiclassical approximation for the propagation kernel in (j4.13p between distinct 
points X and a;', in a time interval t: 



(x\e^"*\x'' 



'det 



52 i? 



dx dx' 



(4.15) 



Here, i?(x, x'; t) is Hamilton's principal function for the path connecting the points x and x' , in time t. The prefactor is 
the Van Vleck determinant factor, given by variations with respect to the path's endpoints. The next step is to evaluate 
the t integral in (|4.13|) by the stationary phase approximation. Here we note that the exponent is jr[Et-\- R{x, x' \ t)], 
whose stationary point condition, ^ = —E^ fixes t to be the period such that the closed classical path has energy 
equal to E. This means that the stationary point condition implements the Legendre transform from R{x, x'\ t) to the 
action S{x^ x' \ E) associated with the closed classical orbit of energy E. Furthermore, the prefactor from the t integral 

contributes a factor 1/ \J = \J which follows from the other Legendre transform relation ^ = The final 

step towards the Gutzwiller trace formula (|4.14p is to take the coincident point limit, x' x, and to integrate over 
the marked point x, the beginning and end of the orbit. Evaluating this integral by a third and final steepest descents 
approximation forces the closed loop to be periodic. The x integral naturally splits into an integral along the loop, 
and an integral transverse to the loop. The integral along the orbit produces a factor of the period Tp, reflecting 
translation invariance with respect to the parameter of the orbit. In the language of quantum field theory, this is a 
collective coordinate contribution, associated with the zero mode coming from the basic fact from classical mechanics 
that the fluctuation operator has a zero mode: K^^x'j^ = 0. On the other hand, the transverse x integrations give yet 



another determinant prefactor l/^/det 
The remarkable result of Gutzwiller 



SSj is that these three different-looking prefactors [coming successively from 
semiclassical approximations to the quantum mechanical path integral, the t integral, and finally the x trace 
combine into a single determinant prefactor that has a simple geometrical interpretation in phase space: 



all 



dx 



dE^ ,j-S(x,x';E) 



92 s 



dx±dx' 



Tper.s[E] 
^det (1 - J) 



(4.16) 



Here J is the monodromy matrix of the periodic orbit having energy E and period Tp. For a precise definition of J 
see [H, [8^[83, including the Maslov indices, which we have suppressed here for the sake of simplicity. It is sufficient 
here to note that the full prefactor is expressed in terms of a classical invariant of the closed orbit, which characterizes 
the behavior of small deviations from the periodic orbit in phase space. 

These steps we have just sketched for the derivation of Gutzwiller's trace formula are closely analogous to the steps 
required in evaluating the worldline path integral expresion (14. ip for the effective action T[A\. Thus, analogous to 
(|4.14p . we seek the following representation of the imaginary part of the effective action: 



ImF = d^x e-"'^ 

e-Sp(m2) 

^ or^P v/dct(l-Jp) 



X e 



-(-Dl)T 



X) 



(4.17) 
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FIG. 9: A pictorial representation of the geometrical meaning of the monodromy matrix J appearing in expressions (|4.16|l and 
(|4.20p . We consider closed orbit [solid, grey line] in phase space, and displace the starting point slightly in the plane transverse 
to the loop, and then propagate according to the equations of motion, returning after one period to another point {Sq , Sp') in 
the same transverse plane. The monodromy matrix relates the starting and end points: {5q',5p') — J {Sq,5p). 



The effective action T[A] in (|4.ip involves a Euclidean path integral, in four dimensional space rather than a path 
integral in three dimensional nonrelativistic quantum mechanics. There is also an extra factor of 1/T because r[^] is 
a log determinant. Nevertheless, the ideas follow through, and in the end we can combine the three prefactors found 
by following the worldline instanton strategy enumerated above, to form one simple prefactor expressed in terms of 
the monodromy matrix for the periodic orbit, now viewed in phase space [s^] . The successive approximations are first 
for the the quantum mechanical path integral: 



K{x,x';T) := (x| 



x\e 
1 



(4.18) 



det 



dx dx' 



'R{x,x';T) 



where R{x,x';T) is the Hamilton principal function for the classical trajectory from x to x' in four-dimensional 
Euclidean space, in the proper-time interval T. Next, the T integral is evaluated by steepest descents. The critical 
point of the exponential factor arises when ^ = —m?. This has a natural classical interpretation in terms of the 
Legendre transformation between the Hamilton principal function R{x,x' ]T) [expressed in terms of the total time 
elapsed along the trajectory] and the action S{x 
R{x,x';T) = S{x,x';E) - ET 



It follows that If 



'; E) [expressed in terms of the constant energy of the trajectory]: 
T. Thus, the critical point of the T integral 



-E. 



and 



occurs when E — 



so that 



e-™ ^K{x,x';T) 



(2^)2 



det 



/ d^R 
\ dx dx' 



2n 



\ d^R \ 
I 9T2 I 



.-Six 



(4.19) 



The final step is the coincident limit x 



ImF 



x^''^ and the trace over x'^*^\ yielding 



v/det (1 - J) 



(4.20) 



The main advantage of expressing the computation in this language of the Gutzwiller trace formula is that the 
total prefactor is encapsulated in a sin gle determinant, which moreover has a natural mathematical and geometrical 
meaning in the Euclidean phase space (SSj . In previous work [38, 39;, Jl, [tI, [13 the various prefactor contributions 
have been evaluated separately, and then combined at the end. Thus, the computational strategy of the worldline 
instanton approach can be stated succinctly as follows: 

1. solve the classical equations of motion in four dimensional Euclidean space to find all closed periodic trajectories 
of " energy" E = — m : such a solution is a "worldline instanton" . 



16 



2. evaluate the classical action Sci{E 



ImF is e 



^) on these trajectories. Then the leading exponential contribution to 



3. compute the prefactor from the monodromy matrix J for the dominant trajectory (ies). 

As an illustration of this procedure, consider again the case of a one-dimensional inhomogeneity, for example, the 
case of a time dependent electric field directed in the X3 direction. As noted above, for such a configuration the 
vacuum pair production probability can be computed in several other ways using WKB methods [si, IH, [79j . 
and so it is possible to make a direct comparison. We can choose a Euclidean gauge field Aslx^) = §/(a;x4), 
where £ characterizes the overall magnitude of the associated electric field, oj characterizes the scale of the time 
dependence, and f{ioxi) is some smooth function. For example, for a constant electric field £{t) = £, we have 
f{x) = x; for a sinusoidal electric field £{t) = £ cos{ujt), we have f{x) = sinh(a;); and for a single-pulse electric field 
£{t) = £sech^(wi), we have f{x) — tan(x). Then the classical action on a periodic trajectory of energy E can be 



written [here y 



e£ 



S{E) 



dxA 



e£ 



f{uJ Xi] 



2E 
7£ 



(4.21) 



Evaluated at £' = m?, this is precisely the exponent, 5(7^) m^7r/(e£), appearing in the standard WKB result for 
the pair production rate jssl . [33 . 0, Fzoj . To evaluate the prefactor, we can choose X4 as xy. Then the transverse 
xz direction is in fact an invariant "flat" direction, so we do not need to perform the transverse integration. This 
illustrates the important point that (|4.20p must be interpreted appropriately when there are physical zero modes 



associated with symmetries of the field. Thus, we go back to (I4.19P and observe that 



-1/ 



aW 



Furthermore, 



the other determinant factor in (|4.19p is easily computed (see [73]) using the Gel'fand-Yaglom formula: 



det 



\ dx dx' 

Thus, relative to the constant spatial volume V3, 



y'dW) 



(4.22) 



ImFWLi 



2(47r)2m 



-S(E} 



(4.23) 



Note that (|4.23p agrees precisely with the conventional WKB result in (|3.10p . with the added semiclassical interpre- 
tation of the various terms. 



V. CONCLUSIONS 



I conclude by reiterating that the ELI project has great promise to open up an entirely new non-perturbative regime 
of QED, and of quantum field theories in general. There are many experimental and theoretical challenges ahead, but 
the future appears bright. Theoretically, the biggest challenge in the non-perturbative arena is to develop efficient 
techniques, both analytical and numerical, for computing the effective action and related quantities, in external fields 
that realistically represent the experimental laser configurations. A lot of progress has been made in this direction, but 
new ideas and methods are still needed. Another important issue is that of back-reaction, which is typically ignored in 
the limit where the electric field strength is below the critical field in (|1.4p . But the combination of an ELI laser with 
an accelerated electron beam may produce effective field strengths well above the critical field, in which case back- 
reaction becomes significant. Insights from these nonperturbative QED studies should have implications for related 
questions in gravitational systems of quantum field theory in strongly curved spaces, and also for non-pertubative 
effects in nonabelian gauge theories [89|. 
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